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Abstract - We discuss consequences of a recent observation that the sequence of 
, periodic orbits in a chaotic billiard behaves like a poissonian stochastic process on 

■ small scales. This enables the semiclassical form factor K sc (t) to agree with predic- 

, tions of random matrix theories for other than infinitesimal t in the semiclassical 

limit. 
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: 1 Introduction 

The spectral form factor K(t), defined as the Fourier transform of the spectral autocorrelation 
function, has a central role in semiclassical analysis of chaotic systems. All spectral statistics 
O ■ bilinear in the density, such as spectral rigidity (A3) and number variance (S2) may be expressed 

in terms of the form factor fil, |). The form factor is a convenient tool for semiclassical analysis 
while it can be expressed as a double sum over periodic orbits B. The approximation of the 
form factor thus obtained K sc (t) is called the semiclassical form factor. This paper contains 
some complementary results to Berry's classical paper where he shows that, in the limit 
of the small r, the form factor is K sc (t) — t for system with no time reversal symmetry and 
K sc {t) — It for system with time reversal symmetry. This agrees with predictions from random 
matrix theories || for the Gaussian Unitary Ensemble (GUE) and Gaussian Orthogonal Ensemble 
(GOE) respectively. This result is essentially the only semiclassical support for universality in 
level statistics for chaotic systems obtained so far. However, as we will see, the result is only 
obvious for smaller and smaller r when the semiclassical limit is approached, and we need a 
mechanism to ensure validity up to fairly large r (of the order of unity) for universality to be 
achieved. 

In this paper we are going to describe the sequence of periodic orbit in a statistical language 
and use tools from the thermodynamic formalism of chaotic systems 0]. These tools will be 
worked out in section 2. The main point is assumption B, which states that the sequence of 
periodic orbits may be described as a poissonian process, at least on small scales. In section 3 we 
show that, under this assumption, Berry's result K — (2)r may be recovered without assuming 
conditional convergence of the Gutzwiller formula S. And even more important, this result holds 
for other than infinitesimal r in the semiclassical limit. We do not discuss when assumption B 
breaks down and the corresponding deviation from K = (2)r may start. In section 4 we discuss 
the various assumptions behind our result and how the results depend on them. In particular we 
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review the numerical evidence for assumption B and in section 5 we discuss related work in view 
of our achievements. 



2 Preliminaries 

The key step is the analytical continuation of the Gutzwiller formula for the level density down 
to the real energy axis. We are going to perform this for the Fourier transform of the level density 
under some simplifying assumptions. 



2.1 The Fourier transform of the level density 

Our starting point is Gutzwiller- Voros zeta function j|, |(| which for systems with two degrees of 
freedom reads 

°° ( P iS p (E) \ 

where A p is the expanding eigenvalue of the Jacobian along the flow. To make things as 
transparent as possible we restrict ourselves to billiards, so that the action integral is simply 
S P (E) = l p ■ k(E) where the 'momentum' is k — \/2E. The level density to be discussed is 
measured in K-space and not, as is usual, in S-space. We use units such that m = h = 1. We 
have neglected the Maslov indices, this neglect is discussed in sections 2.2 and 4. 

It is expected that the zeros of the zeta function approximate the quantum eigenvalues. How 
the the semiclassical approximation and the stationary phase approximation behind Zqv affect 
the spectrum is not all clear. However, this is not our main concern here. What is important 
in this paper is that the zeta function may contain additional zeros reflecting its convergence 
properties. 

In this paper we will use the following simplifying assumption: 
Assumption A: 

The zeta function Z(k) is entire and has zeros on the exact quantum positions ±/tj, i ^ and 
one extra zero kq = on the positive imaginary axis lying on the border of convergence. The 

motivation and discussion of this assumption is postponed until sec 4. (The reason behind the 
choice of subscript of h\ji will be obvious in section 2.2.) 

The level density is usually split up into the mean density and an oscillating part d(n) = 
do(K.) +d osc (K). The leading part of do is the Weyl term do ~ Ak/2it where A is the billiard area. 
The Gutzwiller- Voros zeta function above is derived from the Gutzwiller formula for d osc . For 
later purposes we need the Fourier transform of the oscillating part of the level density 

1 i-oo+iC j 

D osc (l) = — e~ iKl — \ogZ( K )d K . (2) 
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We will evaluate it in two ways. First, by inserting the product representation of the zeta function. 
It is thus essential the the constant C in eq (|J) is sufficiently large so that the contour runs in 
the region where the product (|l|) converges (and thus well above all zeros of Z). One can then 
exchange summation and integration and establish the result. 

^c(/) = E^E l Ml"/ii ■ < 3 > 

p n =l 1 P 1 

Exponential divergence of this sum D osc (l) —> exp(ft, 1 / 2 0i a feature to be discussed in section 2.2, 
is directly related to the presence of the zero k = ihi/2 as can be seen from eq. ^) 
Secondly we compute D osc by means of residue calculus 

oo 

D osc {l) = e h °° l + ■ ( 4 ) 
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The reader may wonder why the contribution from the mean distribution do has disappeared in eq 
. The reason is that we have neglected the contribution from the large semi-circle which would 
have yielded delta functions, and even more nasty things, associated with the Fourier transform 
of do. The tilde-sign in D osc indicates that the result carries contribution from the extra zero kq. 
Without this zero we get the Fourier transform of the true level density D osc — Y^jL-oc j^o e ~ lKjl ■ 
We can then establish the identity 

^c(0=£* P £rg^-e^ . (5) 

p n=l ' P ' 

We have all the way assumed that I > 0. 

From now on we will tacitly replace every occurrence of |M™ — I\ with |A p |. The errors thus 
induced are completely negligible in the limits we are going to explore. We must now discuss 
some properties of the periodic orbits sums. 



2.2 Some properties of the set of periodic orbits 

Essentially all dynamical information is encoded in the sequence of the invariants of the primitive 
periodic orbits {l p , A p } (the maslov indices provide some topological information though). In the 
asymptotic limit I — > oo one can establish the following family of sum-rules for chaotic systems 

E^Ehtf^^ ■ (6) 

(From now on we will reside in the asymptotic limit of large I and write equality signs instead 
of arrows.) The result applies after appropriate smearing of the delta functions. The entropy-like 
quantity hp decreases with increasing (3. The special case (3=1 was discussed already in Qj; for 
bound systems one have hi = 0. ho is the topological entropy ||. General (3 are discussed in e.g. 
§0,0. Usually there are some restrictions on (3. For e.g. the Sinai billiard eq (^|) is only valid 
provided that — 1 < j3 < 1 p0| . In our considerations it suffices if eq. (||) is valid for 1/2 < j3 < 1. 
We expect this property to hold for any reasonable system. 

It is very useful to describe the set of periodic orbits in a statistical language. For large / the 
repetitions of shorter orbits are overwhelmed by the number of primitive orbits so that we may 
neglect the sum over n above. Let us call the density of of prime orbits (j>(l) = ^2 p 6(1 — l p ). From 
eq. (H) we see that the mean value of this density is 

< 4>{l) >= e hal /l . (7) 
For general (3 we express the sum rules in terms of the averages < |A p |~^ > 

I <(/>>< \A\-P >= e h " 1 < lAI^ >= e hlil . (8) 
For instance we get the result which will be of use later 

< |A|- X >=e-' l ° l . (9) 

So much for the large scale structure of the sequence {/ p ,A p }, i.e. large smearing widths in eq 
(||). What about the small scale structure? 

Let us order this sequence according to increasing l p and consider the ordered sequence {li, A^} 
where the integer i denote the position in the sequence. Then rescale the length variable according 
to 

rk 

< <j>{l') > dl' , (10) 
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so that the mean spacing < li — &i-x > is unity. We will now make our main assumption. 
Assumption B: 
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1 / The sequence £j is given by a poissonian process with unit intensity. 

2/ The corresponding stabilities A$ may be considered as mutually independent stochastic variables 
We discuss the evidence for this assumption in sec 4. 

We can now reformulate D osc in a purely statistical language 

D osc (l) = I 1 0(1) 1 -<0><_ L>| . (11) 

The introduction of phase indices (Maslov indices and symmetry indices |]l3| ) will generally 
move down the leading zero hi/ 2 [|llj and there is a possibility that it might even cross the real 
K-axis making the Gutzwiller sum conditionally convergent |fl2f . It is nontrivial to deduce if this 
really takes place for a given system. Many estimations of the position of the entropy barrier, 
with or without Maslovs, in the literature assumes uniform hyperbolicity of the system and are 
thus invalid for generic systems. 

3 The form factor 

The spectral form factor is defined as the Fourier transform of the spectral autocorrelation func- 
tion 

1 f°° 

K=-r dee~ ld d osc ( K + e/2)d osc ( K -e/2) . (12) 

"0 J-oo 

It is usually regarded as a function of the dimensionless length variable r = l/(2ndo) = I /(An) 
with k as a parameter. The suggested universal behaviour of K(t) should arise in the semiclassical 
limit k — ► oo. To be meaningful the form factor needs some averaging which we will apply first at 
the end. The form factor can be expressed in terms of the Fourier transform of d osc (K) according 
to 

1 f°° f°° I' + 1" 

K=— dl 1 dl"6(l - ——)cos(k(1' - l"))D osc (l')D osc (l") . (13) 
2irda J_ oc J_ OQ 2 

The derivation is straightforward, one has to use the fact that d oac (n), and thus D osc (l), are real 
and even. 

Assuming conditional convergence 

Let us now follow Berry's arguments a little longer Q|. If the Gutzwiller formula is conditionally 
convergent we can insert D osc directly instead of D osc into (|l3|). 

K - - 4 / / < "" c ™ w, ' - l " m - ^ )( ? w' ■ (14) 

Keeping k constant and letting I — > be small the cosine will wash away the non diagonal terms 
(assuming no systematic degeneracies among the li due to e.g. time reversal symmetry) and we 
get 

p 6(1 -k) 

K - = ^r %^KT ) ■ (15) 

Using the sum rules in sec 2.2 we get 

^ = ^ r yexp(-M) = ^r=r . (16) 

This result gave rise to some enthusiasm since it agrees with predictions of random matrix 
theories. But one may ask two questions. First, is eq ( |l6| ) true even if the Gutzillwer sum is not 
conditionally convergent? 

Secondly, we note that the number of periodic orbits contained in one period of the cosine in 
eq ( [HI) is ^< > /k ~ exp(hoAKT)/(AK 2 r). To make this estimate we have assumed that the 
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smearing width is Al «< <j> > _1 which is the smallest conceivable choice. In the semiclassical 
limit (k — > oo) the argument leading to eq ( |l6| ) is brutally violated for other than infinitesimal r. 
If the predictions of random matrix theories are correct one expects eq ( filf ) to hold up to t ~ 1. 
Now to the second question. Is the result K — r correct for finite r, and if the answer is yes, 
why? 



The generic case 

The key lies in the stochastic nature of the periodic orbits as formulated in assumption B. First 
we insert the general formula for D osc (l) 

<s< < * x jtn »C < ♦ >< Tbr » ■ ( ' 

We clearly see that we are dealing with correlation in the sequence {k,Ai} and, according to 
assumption B, there arc no correlations at all. Therefore only the diagonal terms will contribute. 
In order to avoid delta functions in the resulting form factor we smear it 
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K sc = ^J K(l')dl' , (18) 
which now may be expressed as 

K sc = -^— f i+A <fl'((r^J-<^X^>)l . (19) 



The integral in this expression is just the variance of the sum of 1/ y|A| in a window of a poissonian 
process, 

I 2 



K 



V A (£^=) . (20) 
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The calculation of this variance is an elementary exercise in probability theory (see Appendix) 

2 

and the result is Va = A A < l/-v/|A| >. The intensity A equals the mean density of prime orbits 
A =< (f)(1) >= exp(hol)/l and according to the sum rules in sec 2.2 we have 

. 2 

< 1/v |A| >=< 1/A >= exp(-/i Z). We thus get our final result 

which is the same as for conditionally convergent systems but the result now holds for other than 
infinitesimal r in the semiclassical limit. 

It is straightforward to generalize to the systematic degeneracies of periodic orbit exhibited 
by time reversible systems, and we will not discuss it. 



4 Motivations for our basic assumptions 

Our result relies on a series of assumptions and approximations. Some of them may be removed 
or modified without altering the result. In this section we motivate our assumptions and discuss 
the extent to which the results depend on them. 

First assumption A. The presence of a leading zero kq — ih\/2 such that h\/2 > (no 
Maslovs) has already been discussed. In the general case it is reasonable to assume the presence 
of several zeros not associated with any quantum state. Their contribution is naturally included 
into < 4> >< | vV | 1 / 2 (Z ) > giving rise to oscillatory and exponentially decreasing corrections. 
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The assumption that the semiclassical zeros equals the quantum eigenvalues is not crucial for 
our results. It was mainly introduced for computational and notational convenience. However, 
the expected failure of the semiclassical eigenvalues to be real will have consequences for the large 
r behaviour of the semiclassical form factor ^0|, see section 5! 

There is one example for which assumption A is exactly fulfilled and that is compact billiards 
on surfaces of constant negative curvature. These systems are special having zeros on the exact 
quantum positions. But there is also a zero on the imaginary k axis, right on the border of 
convergence |l4) . Indeed there is a zero on the border of convergence of each j-factor in the zeta 
function (0). 

Now to assumption B. In ref fl5| the authors pursued the original idea to consider the spectrum 
of lengths lj of the prime cycles and do level statistics a la quantum chaos. Their system was 
a touching three-disk billiard. The first step is to unfold the spectrum, cf section 2.2, yielding 
the sequence lj. Then they studied the level spacing distribution, which was found to be an 
exponential to very high accuracy, and spectral rigidity, which was found to agree with A3 (L) = 
-L/15. This is consistent with the sequence lj being given by a poissonian process and motivates 
assumption B. The three-disk billiard, as well as any generic Euclidean chaotic billiard, has almost 
certainly an infinite symbolic dynamics. It would be nice to know if the greater regularity among 
the cycles for a finite symbolic dynamics still exhibits this kind of randomness. 

We do not attempt to explain this random feature, we only offer the following hand- waving 
argument. Neighbours in the sequence {lj} need not be close in phase space and there is therefore 
no reason for correlations. In a small proximity to a given length I there may be many periodic 
orbits. The sequence of periodic orbit could thus, perhaps, be viewed upon as the random 
superposition of many sequences, each sequence corresponding to a topologically distinct family 
of periodic orbits. This would give rise to the poissonian nature. 

Assumption B2, concerning the the stabilities, is in the same spirit as Bl. If the lengths are 
uncorrelated so should the eigenvalues be. 

The restriction to billiards is mostly for convenience. We find it highly unlikely that a smooth 
potential would exhibit complete chaos. Our calculations would be easily modified for a chaotic 
smooth potential, if existing, which is homogeneous; the 'momentum' variable k being some other 
power of energy E. If one chooses to consider non homogeneous potentials one must perform the 
Fourier transform with respect to Ti instead of k We won't speculate about this case since it would 
take us to far from the case where assumption B has been verified. 



5 Discussion 

It is interesting to note that the semiclassical form factor K sc (r) in the region < r < 1 depend 
on both the very small scale structure and the very coarse structure on the sequence of periodic 
orbits. In this paper we have focused on the deep asymptotic limit k — > 00. For finite energies 
one may have to consider pre-asymptotic behaviour and power law corrections of the periodic 
orbit sum rules involved. This is discussed in r efs H0. In these papers we did not correct for 
the divergence of the trace formula but this procedure is readily justified from the present results. 
This pre-asymptotic behaviour extend considerably the non universal regime in spectral statistics 
derivable form the form factor as compared to the regime discussed by Berry. Such non-universal 
regimes have been observed in several numerical experiments jl7[ |l8|, [if], |ll[ . In ref [|ll| we also 
discuss the role of marginally stable orbits which plays a major role for moderate energies. 

The semiclassical status of the large r limit is much more obscure. It is clear that the 
quantum form factor approaches unity in this limit K — > 1 provided that there are no systematic 
degeneracies in the spectrum. In pQ] Keating demonstrates that, since we cannot expect the 
trace formula to produce poles exactly on the real axis, the semiclassical form factor should 
diverge exponentially. If the imaginary part of the poles is much less than the mean spacing, as 
is indicated by |2]J| , this exponential take off should not occur until fairly large r and one can 
still hope for saturation of the semiclassical form-factor. 

Some evidence of saturation is presented in ref. |2^] for the hyperbola billiard and other 
systems. 



6 



The exponential collapse of K sc (t) reported in |g3| appears to be due to neglect of the 
extra zero(s) and illustrates the hazard of inserting diverging series into the form-factor. This is 
particularly dangerous since the expression thus obtained need not be divergent. 

Whatever happens to the form factor we can of course not expect assumption B to hold 
throughout the spectrum and a deviation form K sc (t) = (2)r should be expected. If the expected 
saturation indeed takes place it is, of course, highly desirably to understand its classical origin 
and its manifestation by the periodic orbits and the connection to the fact that the system is 
bound. 

The reader may think that our assumption B contradicts the concept of action repulsion 
discussed by Argaman et.al. p2fl . However, if the periodic orbit correlations proposed in ref. jp2| 
really occur this repulsion is not an effect acting among neighbours in the sequence {U} of cycles 
but over vast distances, the name action repulsion is thus not very appropriate. 

Much of the present work on the semiclassical trace formula is concerned with taming the 
diverging trace formula and the computation of corrections in order to obtain accurate results for 
the bottom part of the spectrum. Many of these corrections disappear in the semiclassical limit. 
In this paper we took the opposite point of view. We tried to relate asymptotic (=semiclassical) 
properties of the spectrum to the asymptotic behaviour of the periodic orbits. As both are con- 
veniently expressed in a statistical language this approach aims at unveiling the duality between 
periodic orbit statistics and level statistics. 
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Appendix 

Ake Svensson owns a small shop in the old town of Stockholm. Customers arrive at the shop 
according to a poissonian process with intensity A. The amount of money paid by one customer is 
considered as a stochastic variable x with probability distribution f(x). x belonging to different 
customers are mutually independent 

During time T the amount of cash received by Ake is Xt- What is the variance of Xt ? - 

The distribution of arrivals during time T in a poissonian process is 

Pn=e -AT(Arr (22) 

n! 

so the distribution of the variable X is 

oo 

F{X) = Y,Pnt n {X) , (23) 

71=0 

where /*" is the n-fold convolution of /. The mean and variance in such a convolution are additive 
so we have 

< x n >= n < x > (24) 
< x 2 n >= n < x 2 > +n(n - 1) < x > 2 . (25) 

A short calculation now yields the mean and variance of Xt 

< X T >= \T-<x> (26) 
V(X T ) =<X T >- <X T > 2 = XT- < x 2 > . (27) 
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